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We describe a technique coupling standard rheology and ultrasonic imaging with promising applications to 
characterization of soft materials under shear. Plane wave imaging using an ultrafast scanner allows to follow 
the local dynamics of fluids sheared between two concentric cylinders with frame rates as high as 10,000 
images per second, while simultaneously monitoring the shear rate, shear stress, and viscosity as a function 
of time. The capacities of this "rheo-ultrasound" instrument arc illustrated on two examples: (i) the classical 
case of the Taylor-Couette instability in a simple viscous fluid and (ii) the unstable shear-banded flow of a 
non-Newtonian wormlike micellar solution. 
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I. INTRODUCTION 

Soft matter is a flourishing field that involves an ever 
increasing number of characterization tools. Among 
them, techniques that allow one to follow the dynam- 
ics of complex fluids under shear are required to bet- 
ter understand the coupling between the shear flow and 
the material microstructure, i.e. its organization at the 
supramolecular scaled. Indeed, non-Hookean and/or non- 
Newtonian features in soft materials typically arise from 
the modification of the microstructure due to deforma- 
tion and flow. In turn, microstructural changes lead 
to modifications of the fluid viscosity and hence of the 
flow field. For instance, the alignment and elongation 
of polymer molecules by shear are well known to lead 
to a reduction of the viscosity, known as shear thin- 
ning, in polymeric solutions. In some cases, due to this 
feedback mechanism between flow and microstructure, a 
homogeneous shear flow may become unstable and give 
way to a heterogeneous flow constituted of several bands 
of different local viscosities. Such a phenomenon, re- 
ferred to as shear banding in the literature, turns out 
to be widespread among very different classes of materi- 
als, ranging from colloidal gels, emulsions, and pastes to 
self- assembled surfactant and polymer systems^r^. Yet, 
shear banding is just one among several examples of com- 
plex flows found in soft materials, which may also display 
fractures in the bulk or apparent slip at the bounding 
walls^"— . The complete understanding of the physics of 
such dynamical and heterogeneous phenomena calls for 
space- and time-resolved investigations of the flow prop- 
erties in well controlled geometries. 

Rheometry is the most widespread tool to characterize 
the flow properties of a soft material in well-controlled 
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geometriesi^ii^. Shear is generally applied by a rotating 
tool attached to a measuring device involving an air or a 
magnetic bearing. The most common shear cell geome- 
tries are (i) the cone-and-plate, where the sample is con- 
fined between a rotating cone and a fixed plate, (ii) the 
plate-and-plate, where the sample is confined between 
a rotating plate and a fixed plate, and (iii) the Taylor- 
Couette geometry, where the sample is confined between 
a fixed outer cylindrical cup and an inner rotating cylin- 
drical bob. Rheological data captured by a rheometer 
include the shear strain 7, the shear rate 7 = d7/di and 
the shear stress a experienced by the material, which are 
respectively computed from the position of the rotating 
tool, its rotation speed and the torque exerted on the 
tool. Standard rheological protocols consist in applying 
either a constant 7 (cr resp.) and to measure the cor- 
responding a (7 resp.), from which the apparent shear 
viscosity rj = a/j is deduced as a function of time, or to 
impose an oscillating 7 or a, from which the viscoelas- 
tic moduli G' and G", known as the elastic and viscous 
moduli, are respectively computed as the in-phase and 
quadrature-phase component of the measured response. 
For all cell geometries the proportionality factors that 
link the rotation speed and the torque to the shear rate 
and to the shear stress are calculated based on various 
assumptions on the flow and on the material. In par- 
ticular, both the flow and the material are assumed to 
remain homogeneous. The reader is referred to standard 
rheology textbooks for more detail o^^'^'^ . Our main point 
is that rheological measurements are essentially blind to 
local microstructmral and/or flow heterogeneities since it 
provides only spatially averaged observables. Therefore 
complementary tools are needed to investigate complex 
flows such as the shear banding or fracture flows men- 
tioned above. 

In the last two decades a number of combined tech- 
niques have been introduced to go beyond standard 
rheology and get simultaneous local information on 
the microstructure and on the flow field. These have 
been recently reviewed in^. Let us mention local 
structural characterization under shear through light^^. 
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neutroni^ and X-ray^ scattering or through birefrin- 
gence measurement a^^'^^ , as well as local velocity mea- 
surements through optical particle tracking^Sr^^^, ultra- 
sonic velocimetry2^, or magnetic resonance imaging^!. 
Among them the ultrasonic speckle velocimetry tech- 
nique (USV) developed by one of us led to a number of 
important findings on shear banding flows thanks to its 
fast temporal rcsolutio n^^i^^ . The aim of the present pa- 
per is to extend our previous one-dimensional technique^^ 
to a combination of rheometry and ultrafast ultrasonic 
imaging leading to a two-dimensional time-resolved char- 
acterization of the shear flow of soft materials. In Sec. [TT] 
we describe the experimental setup used to acquire simul- 
taneously global Theological data and ultrasonic images 
of the sheared material with emphasis on the technical 
features of our ultrasonic technique and on the process- 
ing of the acoustic data. In Scc. lIIII the instrument is first 
calibrated at low shear rates on a Newtonian fluid, a di- 
lute suspension of hollow glass spheres. We further show 
that it can be used to image the toroidal vortices that 
appear above the critical shear rate corresponding to the 
onset of the classical Taylor- Couette instability. Finally 
Sec. IIVI illustrates the capabilities of the technique in the 
case of a semidilute surfactant solution known to present 
both shear-banding and viscoelastic instabilities^!. 



II. EXPERIMENTAL SET UP 

Our apparatus consists of a custom-made ultrasonic 
scanner coupled to a commercial rheometer. It is derived 
from our previous one-dimensional ultrasonic velocime- 
try technique^^. Figure [Ija) shows the general design of 
the various instruments involved in our technique. Be- 
low wc briefly discuss the rheometer and the rheological 
cell, which are both standard in the field of complex fluid 
characterization. Then we focus in detail on the ultra- 
sonic imaging system. 



A. Rheological measurements 

The rheometer is a state-of-the-art stress-imposed 
rheometer with a magnetic air bearing (TA Instruments 
ARG2). It is equipped with a custom Taylor- Couette 
cell of height 50 mm. The radius of the inner rotating 
bob made of PEEK is Ri = 48 mm while the fixed cup 
made of PMMA has a radius R2 = 50 mm, yielding a 
gap e = R2 — Ri = 2 mm. In a cylindrical geometry the 
shear stress a is not perfectly constant across the gap: it 
rather decreases as l/{Ri+r)'^ , where r is the distance to 
the inner bob. In our case this leads to a decrease of cr by 
Sa/a c^2e/Ri ~ 8 % from the (inner) bob to the (outer) 
cup, which can be reasonably neglected so that the stress 
field can be taken as quasi-homogeneous. Also note that 
the bottom of the inner bob is machined as a cone of an- 
gle 2.3° and truncated at 50 ^m from its tip, so that the 
shear rate below the rotating bob matches the average 




FIG. 1. Sketch and picture of the experimental setup, 
(a) Three-dimensional general view, (b) Top view of the gap 
of the Couette cell together with the path of the acoustic beam 
and the various axes and angles defined in the text, (c) Pic- 
ture showing the water tank with the ultrasonic transducer 
facing a smooth, transparent Taylor-Couette geometry with 
dimensions {Ri = 23 mm and R2 = 25 mm) smaller than 
that used in the text. 



shear rate across the gap. This concentric cylinder cell 
terminated by a cone-and-plate type of geometry at the 
bottom is referred to as a Mooney-Couette cell in the rhe- 
ology literaturei^ii^. Moreover it is well-known that the 
physical-chemical properties of the cell walls may affect 
drastically both global rheological properties and local 
flow bchaviors^ i 1 1 1^8,29 ^ particular wall roughness can 
be used to avoid slippage of the sample on the rotating 
bob and/or on the fixed cup. Here a "smooth" (polished) 
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bob is used together with a "rough" (sand-blasted) cup. 

The cell is immersed in a large water tank (total vol- 
ume of 2.2 L) connected to a water bath (Huber Min- 
istat 230-CC-NR) that keeps the temperature constant 
to within 0.1°C [see Fig. [lja,c)]. This water tank is also 
used to provide acoustic coupling with a fair impedance 
matching between the sample contained in the cell and 
the ultrasonic transducer array (sec below) . The temper- 
ature at which the present experiments are conducted is 
20.2° C everywhere except in Sect. IIVI Finally the thick- 
ness of the cup is 5 mm everywhere except for a rect- 
angular region milled to provide a minimal thickness of 
1 mm and thus prevent a too strong attenuation of the 
ultrasound when travelling across the cup. 

When illustrating the technique in Sees. IIIII and IIVI we 
shall only use the rheometer in the "shear startup" mode 
whereby a given steady shear rate is imposed starting 
from rest at time i = and the subsequent stress response 
a(t) is monitored. However it is clear that the technique 
can be adapted to steady shear from any initial shear rate 
other than zero, to "step strain" experiments where a 
given strain is imposed at t = 0, to "creep and recovery" 
experiments where the shear stress is controlled rather 
than the shear rate, or even to oscillatory shear experi- 
ments. In all the experiments discussed below the target 
shear rate is reached by our stress-imposed rheometer us- 
ing a feedback loop on the bob rotation speed. Note that, 
although the characteristic time of the feedback loop is of 
the order of 10 ms, it may take a much longer time before 
the target shear rate is actually reached by the rheome- 
ter due to the coupling between the geometry inertia and 
the fluid viscoelasticity. With the Couette geometry de- 
scribed above, we found this time to be of the order of 
0.5 s in water [see Figs. [Tljb) andlS^b)] and about 3 s in 
viscoclastic wormlike micelles [see Figs. [SJ^b)]. 



B. Ultrasonic scanner and probe 

Our ultrasonic scanner is built upon a custom-made 
multichannel electronics designed for phased-array appli- 
cations by Lecocur Electronique. It consists of 128 inde- 
pendent channels, each made of a transmitter and a re- 
ceiver. The transmitter can be either a "spike" transmit- 
ter that emits short pulses of tunable voltage (from 10 to 
230 V) and width (from 25 ns to 3 /iS, fall time less than 
10 ns) and whose sequence in time is programmable with 
a maximum pulse repetition frequency (/prf) of 20 kHz, 
or a fully programmable analog transmitter that gener- 
ates power-amplified arbitrary waveforms with a band- 
width of 10 MHz and a maximum peak-to-peak ampli- 
tude of 100 V. In the latter case a 4 Mb memory is used 
to store the waveform which is sampled at 80 MHz and 
converted by a 12-bit linear digital-to-analog converter. 
Analog transmitters allow for transmitting high acoustic 
power at central frequencies up to 10 MHz while pulsed 
emission is used whenever a larger emitting bandwidth is 
needed, e.g. for ultrasonic imaging at frequencies larger 



than 10 MHz. In the present paper only pulsed emis- 
sion is used to drive a high-frequency transducer array 
as described below. 

The receiver is composed of a 12-bit linear analog- 
to-digital converter with a sampling frequency /<, = 
160 MHz and a bandwidth of 1 to 30 MHz, an ultra-low 
noise amplifier with a maximum gain of 80 dB, and a 
SDRAM memory of 32 Msamples. The input impedance 
is 50 f2. A real-time processor board is also featured, 
which allows for a real-time processing of ultrasonic data 
(not used in the present work) and for faster data transfer 
to the computer. The ultrasonic scanner is programmed 
under Matlab and the acquired data are stored on the 
PC hard drive for post-processing as described below in 

sec.iiro 

The ultrasonic probe used in the present work is a 
custom-made linear array of 128 piezoelectric transduc- 
ers designed by Imasonic. The transducers are 200 ^m 
wide and spaced by 50 /im, resulting in a total active 
length of 32 mm. They work at a central frequency 
/ = 15 MHz with a bandwidth at -6 dB of about 8 MHz. 
The wavelength in water is thus A = c// ~ 100 /im, 
with c = 1480 ms^^ the sound speed in water at 20°Ci^ 
In the elevation direction, i.e. in the direction perpen- 
dicular to the imaged plane, the transducers have an 
aperture D = 10 mm with a cylindrical shape of radius 
30 mm (see sketches in Fig. [Ij . This leads to a focus- 
ing in the elevation direction at a distance = 30 mm 
from the probe and over a focal spot whose length If 
and width Wf at -6 dB can be estimated respectively as 
If ~ 8\F'^/D^ ~ 7 mm and w/ ~ \F/D ~ 300 fiin. 

In the following the region of interest, namely the gap 
of the Couette cell, will be centered in the focal spot 
and will be thus considered as a thin rectangular slice of 
height 32 mm, length e = 2 mm, and thickness 300 /im. 
In order to detect a non-zero displacement of the sheared 
material, the transducer array is set at an angle 00 ~ 5° 
relative to the normal of the outer cup. After refrac- 
tion in the cup, the actual angle of incidence within 
the sample, noted 0, is close to but different from 0o 
as shown in Fig. [Ijb). This parameter will be precisely 
determined through the calibration procedure described 
in Sect. Ilira 



C. Ultrasonic data processing 

1. Plane wave imaging 

Standard ultrasonic imaging schemes use a series of 
emissions focused at different locations in the sample 
and from which a single image is reconstructed. Doppler 
imaging can then be performed to measure one compo- 
nent of the flow velocitj^. In more recent ultrasonic 
velocimetry approaches successive images are stored and 
cross-correlated as in optical Particle Imaging Velocime- 
try (PIV) to infer the two-dimensional displacement field 
between two images. In general this "EchoPIV" tech- 
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FIG. 2. (a) Raw speckle signal Si{t, z) recorded after a single plane wave emission as a function of time t and transducer position 
z along the array, (b) Corrected speckle signal Si{t,z) after removal of the fixed echoes in Si{t,z) (see text). The dashed lines 
at t = 38.7 and 41.5 fis indicate the limits of the gap as inferred from the calibration procedure described in Sect. IIII"^ The 
signals are coded using linear color levels. Experiment performed in a Newtonian suspension of hollow glass spheres at 1 wt. % 
in water and sheared at 7 = 10 s~^ (see also supplementary movie 1)^. 



nique is based on a standard scanner and yields a max- 
imum frame rate of about 200 Hz which is enough for 
biomedical imaging and diagnosis^ but too slow to fol- 
low transient phenomena or fast dynamics in soft matter. 
High-resolution imaging systems with working frequen- 
cies as high as 30 MHz now allow for "micro EchoPIV" : 
Poiseuille flows in stenotic vessel-mimic phantoms of di- 
ameter 0.6 mm and blood flow in small animals have been 
imaged with a spatial resolution of 60 fim and frame rates 
up to 100 Hz^. 

Our imaging technique relies on the use of "plane wave 
imaging" as introduced by Sandrin et in order to 
increase the temporal resolution of standard ultrasonic 
imaging by a factor of 10 to 100. A single emission of a 
plane wave is generated by firing all the transducers si- 
multaneously. The signal sent to the transducers is con- 
stituted of a single pulse of amplitude 100 V. The echoes 
backscattered by the sample arc amplified with a gain of 
60 dB, recorded, and stored on the on-board memories of 
the various electronic channels. The pulsed plane wave 
emission is repeated N times with a repetition frequency 
/pRF, which constitutes one "sequence." Identical se- 
quences can be repeated up to iVseq = 512 times. The 
number of pulses TV per sequence is limited to 8,192 emis- 
sions and by the total memory of 32 Msamples available 
per channel. Once the full series of imaging sequences 
is completed, the ultrasonic data are downloaded to the 
PC for post-processing. 

Such plane wave imaging has been used to follow shear 
waves into the human body and infer useful informa- 
tion on local elastic properties, a technique known as 
"transient elastography," which was recently taken to 



real-time thanks to graphical processing units^. The 
pulse repetition frequency /prf, hence the frame rate, 
is only limited by the total time-of-flight from the ultra- 
sonic probe to the inner bob and back to the probe. In 
the present experimental setup the total travel distance 
is about 60 mm, which corresponds to a time-of-flight of 
about 40 /iS and to a maximum frame rate of 25 kHz. 

In our case, the aim is to follow the deformation and 
flow of soft materials. The backscattered ultrasonic sig- 
nal either arises from the fluid microstructurc itself or 
is obtained artificially by seeding the sample with mi- 
crospheres that play the role of acoustic contrast agent. 
In both cases we call this signal an "ultrasonic speckle" 
and the corresponding velocimetry technique "ultrasonic 
speckle velocimetry" (USV) . Note that a seeding concen- 
tration of a few volume percent generally provides a good 
signal-to-noise ratio and does not constitute any signif- 
icant limitation of USV when compared to EchoPIV—. 
Still plane wave imaging leads to a strong loss of lateral 
resolution, which is minimized in the following by a re- 
ception beam-forming, so that speckle tracking in USV 
is usually performed only on the axial component of the 
velocity. 



2. From raw speckle data to ultrasonic images 

a. Raw speckle signal. - The ultrasonic speckle 
backscattered from the sample to the transducer array 
results from the interference of all scattering events of 
the incident plane wave by the scatterers across the in- 
sonifled slice (32 mm x 2 mm x 300 /im). If one can 
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neglect echoes that are multiply scattered along their 
propagation, there is a direct correspondence between 
the time-of-flight of the ultrasound and the position of 
the scatterer— . Note that, as in any imaging technique, 
multiple scattering constitutes a strong limitation of the 
present instrument, especially in the case of concentrated 
suspensions, for which a specific analysis of the multiply 
scattered sound should be undertaken^. In the opposite 
case of a sample that is transparent to ultrasound, seed- 
ing the material with acoustic contrast agents allows for 
a good control of the scattered signal. 

In any case the signal recorded on a given reception 
channel is the sum of the echoes from scatterers at vari- 
ous locations in the imaging plane. An example of such 
a raw ultrasonic "speckle" signal recorded in a dilute 
aqueous suspension of polydisperse hollow glass spheres 
(Potters Sphericel, mean diameter 6 /xm, mean density 
1.1 g.cm"'^) is shown as a function of both space and 
time in Fig. [2][a). Here, the time t corresponds to the 
time at which the pressure signal is sampled by the re- 
ceiver, with the origin t = being the time at which the 
incident pulse is sent. In other words t is the ultrasonic 
time-of-fiight from and back to the transducer array. The 
vertical coordinate z denotes the position along the trans- 
ducer array, the origin being taken at the upper end of 
the array, which is about 6 mm from the top of the Cou- 
ette cell [see Fig. (Ha)]. This two-dimensional represen- 
tation of the radio-frequency (RF) backscattered signals 
s{t,z), also referred to as a "B-scan" in the ultrasound 
literature, clearly shows the echoes of the incident plane 
wave on the outer cup and on the inner bob as vertical 
lines at t ~ 38.7 and 41.5 /is respectively. In between 
the bob and the cup, the RF signal is constituted of a 
superposition of spherical wavefronts that correspond to 
the backscattered echoes from the glass spheres within 
the gap. 

b. Procedure for removing fixed echoes. - A close 
investigation of a series of N successive RF signals 
{si(i, jv corresponding to N successive incident 

plane waves reveals that the spherical wavefronts from 
the glass spheres, which move along with the scatterers 
as they are carried away by the shear flow, are mixed 
with other wavefronts that remain fixed under shear (see 
supplementary movie 1 with N = 200)^. These fixed 
echoes are due to waves diffracted by the edges of the 
transducer array which get reflected on the various sur- 
faces surrounding the Couettc cell (free surface, bottom 
of the water tank, Couette cell, etc.). Most of these re- 
flections occur before the main plane wave enters the gap 
of the Couette cell {t ~ 38.7 /is) but the corresponding 
cylindrical echoes may extend over the whole temporal 
window of interest, as seen in Fig. ^a) (see also sup- 
plementary movie 1)^. Moreover the roughness of the 
outer cup as well as small defects in the inner bob may 
also generate fixed spurious echoes. 

As in biomedical imaging where fixed echoes due to 
artery walls may lead to large artifacts in blood speed 
detection^^i^, we get rid of these spurious echoes by sub- 
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FIG. 3. (a) Beam-formed image Si{y,z) computed from the 
corrected speckle signal of Fig. [IJb) and shown as a function 
of y, the distance to tlie transducer array, and z, the position 
along the array. Si is normalized by its maximum value and 
coded in linear color levels, (b) Two successive beam-formed 
speckle signals Si(y,z) (in black) and Si+i{y,z) (in red) for 
a given position z — 14.9 mm along the transducer array. 
Si and Si+i correspond to two different plane wave emissions 
separated by (5t = 2 ms. (c) Enlargement of Si{y, z) (in black) 
and Si+i{y, z) (in red) over a window of width Ay = 2A close 
to the inner bob [indicated as a dashed box in (b)]. This ev- 
idences a noticeable displacement of the speckle to the right 
when going from Si to Si+i. The dashed lines at j/ = 28.7 
and 30.7 mm indicate the limits of the gap as inferred from 
the calibration procedure described in Sect. IIII Al Same ex- 
periment as in Fig. [2] (see also supplementary movie 2)^. 



trading the speckle signal {si{t, z))i averaged over TV suc- 
cessive signals to each individual signal Si{t, z): 



i(t,z) = Si{t,z) ~ {Si{t,z))i 



(1) 



The result obtained from Fig. [2]^a) with a series of = 
200 pulses separated by (5t = 2 ms (i.e. a pulse repeti- 
tion frequency /prf = 500 Hz) is shown in Fig. [21[b) . In 
this case, the Newtonian suspension is sheared at a con- 
stant shear rate 7 = 10 s~^. Shear is started long enough 
before collecting the ultrasonic data so that the fiow is 
laminar and stationary. Figure [2Ub) shows that our pro- 
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FIG. 4. (a) Displacement map 5yi(y, z) computed from two successive beam-formed images Si and Si+i separated hy 5t = 2 ms. 
(b) Displacement map {5yi{y, z))i averaged over 199 correlations between 200 successive images. The dashed lines at y = 28.7 
and 30.7 mm indicate the limits of the gap as inferred from the calibration procedure described in Sect. IIII"^ Same experiment 
as in Fig. [5] 



cedure for removing fixed echoes is very efficient. It re- 
mains so as long as all scatterers across the gap move by 
a significant quantity over the total duration of the pulse 
sequence so that the (moving) backscattered wavefronts 
cancel out in the average. In some cases where the fluid 
may remain fixed in parts of the gap (e.g. in yield stress 
fluids that show shear localization), the average is taken 
on a different series of incident pulses performed under a 
high shear rate, so that the whole sample flows. Since the 
USV algorithm has already been described as a tool for 
flow visualization in other contexts^SilL, we briefly recall 
the main steps of image formation and processing in the 
following and emphasize only the aspects specific to the 
case of cylindrical Couctte flow imaging. 

c. Computation of beam-formed images. - In this 
data processing step based on wave propagation in a 
single scattering regime, an ultrasonic image Si{y,z) is 
formed from a corrected speckle signal Si{t,z) by using 
standard parallel beam forming^^: 

St{y,z) = ^Si{t{y,z,zo),z) , (2) 

where y is the distance from the transducer array along 
the ultrasonic propagation axis (or "axial distance" ) and 

t{y, z, zo) = (3) 

c 

is the ultrasonic time-of-flight from a scatterer located at 
{y, z) in the gap to the transducer located at height zq, 
c denoting the sound speed in water. Note that Eq. (|3]) 
neither accounts for the sound speed in the outer cup 
which is made out of PMMA (cpmma — 2500 m.s~^) nor 
from the fact that the sound speed in the sample may 



differ from that of water. In any case the actual distance 
r in the sample to the inner bob will be inferred from y 
through a calibration step described in Sect. lIIlXl More- 
over the summation in Eq. ([2]) is taken over 30 values of zq 
around z. This fastens the data processing without sig- 
nificantly degrading the image resolution. It can also be 
noted that backscattered echoes are recorded beyond the 
position of the inner bob (i.e. for t > 41.5 fj,s in Fig.[2]and 
for y > 30.7 mm in Fig. [3]). These echoes correspond to 
the scattering of the reflection of the incident plane wave 
on the hoh-^. Therefore only points with y < 30.8 mm 
are considered in Eq. ([3]). The beam- formed image com- 
puted from Fig.[31^b) is shown in Fig.[31Ja) and the line at 
z = 14.9 mm of two successive individual images Si{y, z) 
and Si+i{y,z) are shown in Fig. EJb) and (c) (see also 
supplementary movie 2 for a series of iV = 200 succes- 
sive images separated by 2 ms in a sheared suspension of 
hollow glass spheres)^. 

3. Displacement field 

The speckle displacement along the y direction at point 
(y, z) between two successive pulses i and i -I- 1 is com- 
puted by looking for the maximum of the following cor- 
relation coefficient as a function of J y23,4o . 

y+Ay/2 

C,{y,z,5y)^ ^ S,{y\ z) S,+i{y' + Sy, z) , (4) 

y'=y-Ay/2 

where Ay = 2A ~ 200 /im. An example of correlation 
window of width Ay is shown in Fig. [3{b) and (c) . For 
a given (?/, z), the value of Sy that maximizes Ci{y, z, Sy) 
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FIG. 5. (a) Axial velocity profiles Vy^i{y,z) computed from 
the two successive beam-formed images 5, and Si+i used in 
Fig. |4ja) and shown for z = 7.4 (♦), 15.1 (•), and 22.4 mm 
(■). (b) Axial velocity profiles Vy{y,z) = {vi{y,z))i averaged 
over 199 correlations between 200 successive images. The grey 
line shows the best linear fit of the full data set averaged over 
z. The dashed lines at y — 28.7 and 30.7 mm indicate the 
limits of the gap as inferred from the calibration procedure 
described in Sect. IIII Al Same experiment as in Fig. [2] 



is computed from a second-order polynomial interpola- 
tion of Ci and yields the axial displacement 5yi of the 
speckle between pulses i and i + 1. In practice, to avoid 
any redundancy in the velocity data, such a displace- 
ment Syi{y,z) is computed only for locations within the 
gap such that yk = yo + kAy/3, where yo corresponds to 
the beginning of the gap on the cup side and k is an inte- 
ger. This yields about 30 measurement points across the 
gap of the Couette cell separated by roughly 65 fim. A 
displacement map computed from two successive beam- 
formed images Si and S'i+i is shown in Fig. |31[a). The 
velocity gradient is clearly visible even if the level of noise 
is rather large because only two successive pulses are con- 
sidered. 

Figure Hljb) presents the displacement map obtained 
by averaging Syi(y, z) over the A'^ — 1 = 199 correlations 
corresponding to the sequence of 200 pulses shown in sup- 
plementary movies 1 and 3^. This 0.4 s time-average 
results in a very small level of noise with a typical vari- 
ation of less than 5 % from one channel to another (see 
also inset of Fig. [5]). As expected the flow is laminar and 
homogeneous over the whole height of the Couette cell. 

Axial velocity maps are easily deduced from the dis- 
placement field by Vy^i{y,z) = Syi{y, z)/6t, where 5t = 
1 //pRF is the time interval between two successive pulses. 
By "axial velocity" Vy we mean the projection of the 
local velocity vector v ~ (vryVg^Vz) on the ultrasonic 
propagation axis y [see Fig. [IJb)]. This obviously as- 
sumes that the scatterers follow the fluid velocity field as 
passive tracers. Three individual axial velocity profiles 
Vy^i{y, z) deduced from Fig.|4l[a) are gathered in Fig.[5lja) 
for three different heights along the transducer array, 
while Fig. [Sfb) shows the corresponding time-averages 
Vy{y,z) = {v,{y,z)),. 
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FIG. 6. Tangential velocity profiles v{r) deduced from the 
calibration procedure [Eq. ^ and ©] with yo = 28.66 mm 
and = 5.2° for different applied shear rates: 7 = 5 (•), 
10 (□), 15 (•), and 20 s^^ (T). Data are averaged over 199 
successive correlations and over the vertical direction z. Error 
bars show the standard deviation over z. Grey lines are the 
theoretical predictions for a Newtonian fluid [Eq. (|10[) ]. Inset; 
relative deviation Sv/v computed as the standard deviation 
of v{r,z) over z relative to its mean value v{r). Experiments 
performed in a Newtonian suspension of hollow glass spheres 
at 1 wt. % in water. 



III. CALIBRATION AND TEST IN A NEWTONIAN 
FLUID 

A. Calibration procedure 

In order to recover the velocity field in sheared com- 
plex fluids, it is necessary to calibrate the axial velocity 
data Vy. Such a calibration procedure has already been 
described in detail in Ref*^ and allows one to infer the 
precise value of the incidence angle and the exact posi- 
tion yo of the cup wall. In brief, the calibration procedure 
must be performed in a Newtonian fluid at a low enough 
shear rate to ensure that the flow is laminar and purely 
tangential so that v = {0,vg,0) and Vy simply measures 
the projection of the orthoradial velocity component vg . 
In our geometry and for a dilute aqueous suspension, this 
corresponds to 7 < 40 s~^. Linear fits of the Vy vs y data 
averaged over N ~ 200 pulses and over the vertical di- 
rection z for different values of 7 yield the position yo 
of the cup wall with an uncertainty of about 20 ^m [see 
the fit in Fig. [SJb) for 7 = 10 s~^]. Moreover the angle 
of incidence (j) is estimated within ±0.1° by looking for 
the value of (j> that best matches the expected velocity 
profiles for various applied shear rates^^. 

Once yo and (j) are calibrated, the axial velocity Vy{y, z) 
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is easily converted to a an apparent tangential velocity 
v{r, z), where r denotes the distance to the rotating inner 
bob, through the relationships: 

r = + (y - yo)2 _ 2i?2(y - 2/o)cOS0 - i?i , (5) 

In the small-gap approximation e/Ri <C 1, the previous 
equations reduce to: 



r ~ e - (?/ - yo)sin</), 
v(r, z) ~ —2^ — — 



(7) 
(8) 



Once again it is crucial to emphasize that w(r, z) deduced 
from Eqs. ([5]) and indeed corresponds to the ortho- 
radial component ve of the velocity field v only when 
the fiow is purely tangential. In the general case of 
a three-dimensional fiow with non-zero radial and ver- 
tical components Vr and u^, the USV measurement is 
Vy ~ vg sin (j) + Vr COS (j) in the small-gap approximation. 
Therefore the velocity v defined by Eqs. ([6]) and ([8]) is 
insensitive to the vertical component but may contain a 
significant contribution from the radial component: 



v{r, z) ~ vg{r, z) + 



Vr{r,z) 



taut 



(9) 



Figure |6] presents the result of the calibration proce- 
dure in a suspension of hollow glass spheres at 1 wt. % 
in water. The ultrasonic data are acquired using a sin- 
gle sequence of = 200 pulses with a pulse repetition 
frequency that is tuned depending on the applied shear 
rate 7 according to /prf = 5O7. Whatever the applied 
shear rate, the velocity profiles almost perfectly coincide 
with the theoretical predictions for a Newtonian fiuid: 



v{r) = f'o ( 1 + 



\Rl+r ) 



1 



0-D 



(10) 

where vq is the velocity of the inner bob imposed by the 
rheometer. The last term in Eq. ((TO)) results from the 
small- gap approximation for which vq ~ 76. The only 
significant deviations from the predictions are observed 
close to the outer cup {r ~ 2 mm) where the velocity 
is slightly overestimated. Such an artifact can be at- 
tributed to our procedure for removing fixed echoes which 
prevents us from measuring vanishingly small velocities 
close to the cup. 



B. Shear start-up in a Newtonian fluid 

In order to highlight the time-resolved capabilities of 
our rheo-ultrasonic technique, we turn to measurements 
performed during fiow build-up in a Newtonian suspen- 
sion after shearing at 7 = 20 s~^ is started at time 



t = 0. In the remainder of this paper t shall denote 
the time after start up of shear (rather than the ultra- 
sonic time-of-flight as in Sect. |n|. Here one sequence of 

= 8, 000 pulses is used with /prf = 1 kHz. Figure[7I^a) 
shows a few velocity maps corresponding to averages over 
50 successive pulses every 25 ms (see also supplemen- 
tary movie 3)^. The shear rate and shear stress signals 
recorded simultaneously by the rheometer are also shown 
in Fig. [Zl^b) . Since the fiow clearly remains homogeneous 
along the vertical direction during the whole transient, 
we average the velocity data along z and plot a few ve- 
locity profiles and time series {v{r,z,t))z in Fig. [7{c). 

As expected, the velocity profile is fully developed only 
for t > T, where t = e'^ /v ~ 1 s is the viscous dissipation 
time on the size of the entire gap. At earlier times, a 
boundary layer develops at the rotating inner bob and 
propagates towards the fixed outer cup. Note that the 
velocity time series v{r,z,t) in Fig. [TJc) resemble but 
do not follow exactly the simple error function solution 
for a boundary layer on an infinite flat plate started im- 
pulsively^2il^. The existence of the outer bob, the cur- 
vature and closure of the streamlines, and the complex 
input function for the velocity of the bob all have a non- 
negligible impact on the transient flow. The complexity 
of the input function 7(f) in the inset of Fig.[7{b) is due to 
the fact that we are using an intrinsically stress-imposed 
rheometer working in strain-imposed mode through a 
feedback loop. In particular the fact that the shear stress 
presents a strong undershoot into negative values can be 
attributed to the coupling between the rheometer feed- 
back on the applied torque, the strong inertia of our 
home-made geometry, and the fluid response. Experi- 
ments performed with a smaller, lighter bob would show 
a faster convergence of the applied shear rate to its tar- 
get value with much smaller oscillations. Working with 
a strain-imposed rheometer would allow a better access 
to the actual fluid response to a step- like shear rate func- 
tion. 



C. Taylor-Couette instability in a Newtonian fluid 

We have seen in the previous section that as long as the 
velocity field remains purely orthoradial, the measure- 
ments at different axial positions do not bring any more 
information than one-dimensional velocimetry: the ve- 
locity maps in Fig. [Tl^a) are essentially invariant along z. 
The situation changes radically at higher imposed shear 
rates when this invariance is broken, revealing the full 
potential of our new technique. For a Newtonian fluid 
sheared in between concentric cylinders Taylor showed 
that the purely orthoradial base fiow becomes unstable 
and gives way to a cellular pattern in which the fluid 
travels in helical paths around the cylinders in layers 
of counter rotating vortices-now known as Taylor vor- 
tices^. In the simplest case where only the inner cylinder 
is rotating and the gap is small, the v ortex flow devel- 
ops for Ta > 41, where Ta = y/e/RiRe is the Taylor 
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Oi 




FIG. 7. Start-up of shear in a Newtonian suspension of iiollow glass spheres at 1 wt. % in water and sheared at 7 = 20 s^^ in the 
laminar regime, (a) Velocity maps v{r, z, t) at different times t indicated on the top row. Each map corresponds to an average 
over 50 pulses sent every millisecond (see also supplementary movie 

3)31. (b) Stress response a{t) recorded simultaneously to 
the velocity maps. The symbols indicate the times corresponding to the images shown in (a). Inset: magnification of the stress 
response (j{t) (in black) together with the instantaneous shear rate ■y{t) (in red) imposed by the rheometer. (c) Velocity profiles 
{w(r, z, t))z averaged over the whole height of the transducer array and shown at f = (□), 0.075 (o), O.f 5 (A), 0.225 (V), 0.375 
(o), 0.7 (<), and 1.275 s (>) consistently with (a) and (b). The grey line shows the velocity profile expected for a Newtonian 
fluid in the laminar regime [Eq. (jlOl) ]. Inset: time evolution of {v{r,z,t)}z at r = 0.01, 0.46, 0.98, 1.43, and 1.95 mm from top 
to bottom. 



number (textbooks often use Ta^) and Re = jr is the 
Reynolds number. In the geometry used in our exper- 
iments, the threshold of the instability therefore corre- 
sponds to 7c — 50 s~^, which we indeed observe. 

Figure [D^a) shows consecutive velocity maps for a 
start-up of shear at 7 = 80 s^^ (see also supplemen- 
tary movie 4)^. The first few instants arc still invariant 
along z and resemble the transient flow described in sec- 
tion [TTlBl However, for t > 1.5 s, the symmetry is pro- 
gressively lost with the emergence of well defined undu- 
lations along z. The onset of the vortex flow can also be 
detected on the transient rheological data in Fig. [5{b). 
Indeed, the transition from the purely orthoradial flow 
to the vortex flow brings in an additional flow resistance 
characterized by an increase in the apparent viscosity 
r](t) = <7{t)/'y(t) recorded by the rheometer, which is ev- 
ident on the inset of Fig. MX>) for 1.5 < t < 2.5 s. 

If the flow is observed on longer time scales, the vortex 
flow is observed to be time periodic as can be checked in 



supplementary movie 5^. This behavior corresponds to 
the "wavy vortex flow" regime, the outcome of the second 
instability of the Taylor-Couette system when only the 
inner cylinder is rotating^. For shear rates closer to 
the threshold jc, we also observed steady vortex flows 
(the so-called "Taylor vortex flow" regime, for jc < j 
66 s~^, as expected^). 

Figure [8ja) shows that the main orthoradial flow is de- 
formed by the influence of the secondary vortex flovs4^. 
Typically, at locations along z that correspond respec- 
tively to the outflow/inflow boundary of vortices (strong 
outward/inward radial flow), the main flow is pushed 
outward/inward. Figure [5{c) gives examples of such de- 
formed individual velocity profiles at the locations of an 
outflow boundary (white symbols), an inflow boundary 
(black symbols) and in between (gray symbols). Thus, 
for each wavelength of the undulations on the velocity 
maps, there is a pair of counter-rotating vortices. The 
large field of view of our imaging technique allows us to 
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FIG. 8. Start-up of shear in a Newtonian suspension of hollow glass spheres at 1 wt. % in water and sheared at 7 = 80 s^^ 
in a vortex flow regime, (a) Velocity maps v{r,z,t) at different times t indicated on the top row. Each map corresponds to 
an average over 50 pulses sent every 0.5 ms (see also supplementary movie 4)^. (b) Stress response a{t) (in black) recorded 
simultaneously to the velocity maps together with the instantaneous shear rate 7(f) (in red) imposed by the rheometer. The 
symbols indicate the times corresponding to the images shown in (a). Inset: apparent viscosity r;(f) = a{t)/'y{t). (c) Velocity 
profiles v{r,z,t) at t = 2.375 s and for z = 19.25 (white symbols, outflow boundary), 20.5 (gray symbols, in between outflow 
and inflow), and 21.5 mm (black symbols, inflow boundary). The gray line shows the velocity profile expected for a Newtonian 
fluid in the laminar regime [Eq. (|10|) ]. 



visualize more than 7 wavelengths, yielding a value of 
4.1 ± 0.1 mm, in agreement with the theoretical predic- 
tion Ac — 2e = 4 mm computed by Taylor—. 

To the best of our knowledge a characterization of the 
Taylor-Couette instability solely based on a time resolved 
velocity map of the orthoradial velocity profiles vg(r, z, t) 
at one single location along 6 does not exist in the lit- 
erature. Wereley and Lueptow rather used PIV to mea- 
sure Vr and Vz in a slice of the gap at a given angle 6^. 
Then from the orthoradial component of the momentum 
balance they could reconstruct orthoradial velocity maps 
similar to those measured here. Later, Akonur and Luep- 
tow used PIV on various slices to reconstruct a three 
dimensional velocity map^. Nonetheless, the measure- 
ments at various locations along z were obtained inde- 
pendently at different times and required phase align- 
ment!^. A similar procedure was used by Raguin and 
Gcorgiadis on magnetic resonance velocimctry data of a 
Taylor-Couette-Poiscuillc flow (i.e. with superimposed 
axial flow)^. In contrast to these earlier results, our new 



ultrafast ultrasonic imaging technique provides a way to 
monitor the flow over a large range of vertical positions 
simultaneously and with a temporal resolution that is 
fast enough to capture unsteady flows. 



Note that according to Eq. ^ the data displayed in 
Fig. [Slja) and (c) are not exactly equal to the orthora- 
dial velocity component. The vertical oscillations on the 
velocity maps are enhanced by the contribution of the 
radial velocity to the measurement. Moreover, the veloc- 
ities measured near the inner and outer cylinders devi- 
ate from the expected values more significantly than in 
the purely orthoradial flow regime used for calibration in 
Fig. [51 Once again, this could be an artefact due to the 
procedure used to remove fixed echoes, or it could also be 
a consequence of the secondary flow. Indeed, the vortices 
create localized jets near the walls that can segregate the 
tracers^. 
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IV. ILLUSTRATION IN A WORMLIKE MICELLAR 
SOLUTION 

To illustrate our new imaging technique on a non- 
Newtonian flow of soft matter, we choose an aqueous 
surfactant solution of cetyltrimethylammonium bromide 
(CTAB) at concentration of 0.3 M and sodium nitrate 
(NaNOa) at concentration of 0.34 M. The temperature is 
fixed at T = 28 ± 0.1°C. In this regime of concentration, 
salinity and temperature, the surfactant molecules ag- 
gregate in long, flexible wormlike micelles that entangle 
in a way similar to classical polymers. Such solutions of 
wormlike micelles have been considered as model systems 
for rheological research for more than twenty years^^— . 
They have a Maxwellian behavior in the linear rheology 
regime (i.e. they have a unique viscoelastic relaxation 
time Ty), and they exhibit shear-banding and viscoelas- 
tic instabilities in the non-linear regime^l. 

For a purely viscoelastic fluid sheared in between con- 
centric cylinders, Larson, Shaqfeh and Muller showed 
that the purely orthoradial base flow becomes unsta- 
ble and is replaced by a cellular pattern reminiscent of 
the Taylor vortices but driven by normal stresses rather 
than by inertia^"—. Derivations from simple viscoelastic 
models like the Upper Convected Maxwell model^ and 
experiments on polymer solutions showed that in this 
case the dimcnsionless number controlling the instability 
is a viscoelastic Taylor number S = -y/ e/RiWi, where 
Wi = 7T„ is the Weissenberg number—. Recently, it was 
shown that shear-banding wormlike micelles solutions 
also exhibit a purely viscoelastic instability^!. The base 
flow is not homogeneous but shear-banded, which means 
that the velocity profile in a shear flow presents two dis- 
tinct slopes defining a band of high shear rate (close 
to the rotating inner bob) next to a band of low shear 
rate ji (close to the fixed outer cup)^. The proportion 
a of the high shear rate band increases roughly linearly 
with the global shear rate between 7; and jfj^^i^^n^L^. 
In this case, the instability originates solely in the high 
shear rate band and Taylor-like vortices develop in this 
band for a viscoelastic Taylor number rescaled on the size 
of the high shear rate band S* = y^ae/RiWih, where 
Wih = 'jhTy is the Weissenberg number in the high shear 
rate band^i. 

Figure EJa) shows consecutive velocity maps for a 
start-up of shear at 7 = 46 ± 4 s~^ (see also supplemen- 
tary movie 6)^. In order to follow the flow dynamics 
over about one minute while keeping up with the mem- 
ory constraints of the ultrasonic scanner, the multiple se- 
quence acquisition procedure described in section III C II 
is used with A^scq = 200 sequences of iV = 100 pulses 
at /pRF = 1 kHz. Sequences are separated by 0.5 s. 
Here, the rather poor control of the applied shear rate 
[to within about 10 %, see Fig. EJb)] is due to the cou- 
pling between the rheometer feedback loop, the strong 
sample viscoelasticity, and the large inertia of our home- 
made Couette geometry. 

During the first few seconds of the transient, the flow 



remains invariant along z and shows the development of 
a shear-banded velocity profile, in accordance with the 
rich literature on the subjecli^. More precisely. Fig. [HJc) 
displays a few velocity profiles {v{r, z,t))z averaged along 
z during the nucleation of the shear band for t < 10 s. 
The evolution of these velocity profiles is fully consistent 
with previous works2i^i2i. Note that the establishment 
of the shear-banded velocity profiles is also associated 
with an enhanced wall slip at the moving inner cylinder. 
This fact seems to be independent of the existence of 
secondary flows, as already reported irt^ i^^i^" . 

At later stages (for t > 10 s), the vortex flow emerges 
and deforms the interface between the shear bands. In 
the asymptotic state (for t > 35 s), for each wavelength 
of the undulations on the velocity map, there is a pair 
of counter- rotating vortices, the size of which scales with 
the width of the high shear rate band^l. As shown in 
the inset of Fig. ^h), the emergence of the vortex flow 
generates a viscosity increase, which is already well doc- 
umented^l and is reminiscent of that seen in Fig. [Sl^b). 
As in the Newtonian case described in section UlI CI the 
emergence of the vortex flow is monitored through its 
impact on the main orthoradial flow [see rightmost im- 
ages in Fig. m^a)]. In the inset of Fig. EJc), examples 
of shear-banded velocity proflles deformed by the sec- 
ondary vortex flow are shown at the locations for an 
outflow boundary (white symbols), an inflow boundary 
(black symbols) and in between (gray symbols). When 
compared to the Newtonian case shown in Fig. [HJc) , indi- 
vidual velocity profiles in wormlike micelles are affected 
by the secondary flow only in the highly sheared region 
(for r < 1 mm) since vortices develop only in the high 
shear band. 



V. CONCLUSION AND PERSPECTIVES 

The above results show that our new ultrafast ultra- 
sonic imaging technique is a powerful tool for measur- 
ing spatiotemporal variations of the main flow of sheared 
fluids, simultaneously to traditional rheometry. We have 
successfully tested the technique on a Newtonian fluid 
that undergoes a purely inertial Taylor-Couette instabil- 
ity, and on a non-Newtonian fluid that exhibits shear- 
banding and undergoes a purely elastic Taylor-Couette 
instability. In both cases, our technique can monitor the 
emergence of secondary flows through their impact on 
the main flow, which loses its invariance along the verti- 
cal direction z. Our excellent time resolution gives access 
to the transient establishment of secondary flows and to 
their eventual spatiotemporal dynamics. 

In the past few years the existence of secondary flows 
has been raised as a possible explanation for the often 
complicated signals recorded by rheometers on flows of 
complex fluids, which are usually assumed to be visco- 
mctric. In future work, we plan to use our new imaging 
technique to test for the presence of spatially and tempo- 
rally heterogeneous flows in a variety of complex fluids. 
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FIG. 9. Start-up of shear in a solution of wormlike micelles ([CTAB]=0.3 M, [NaNO3]=0.34 M) seeded with hollow glass 
spheres at 1 wt. % and sheared at 7 = 46 ± 4 s~^ in the shear-banding regime with Taylor-like vortices, (a) Velocity maps 
v{r,z,t) at different times t indicated on the top row. Each map corresponds to an average over 100 pulses sent every 1 ms 
(see also supplementary movie 6)^. Each of the TVsoq = 200 sequences of A*' = 100 pulses is separated from the next one by 
0.5 s. (b) Stress response a{t) (in black) imposed by the rheometer recorded simultaneously to the velocity maps together with 
the instantaneous shear rate 7(t) (in red). The symbols indicate the times corresponding to the images shown in (a). Inset; 
apparent viscosity ri{t) = a{t)/^{t). (c) Velocity profiles {v{r, z,t)) z averaged over the whole height of the transducer array 
and shown at t = 1.8 (□), 2.3 (*), 2.8 {+), and 9.8 s (A). The gray line shows the velocity profile expected for a Newtonian 
fiuid in the laminar regime [Eq. (|10p ]. Inset; velocity profiles v(r,z,t) a.t t — 39.8 s and for z = 10.5 (white symbols, outflow 
boundary), 9.5 (gray symbols, in between outflow and inflow), and 8.5 mm (black symbols, inflow boundary). 



which do not need to be optically transparent but only 
need to scatter ultrasound in the single scattering regime. 
Indeed, our instrument allows us to probe any possible 
flow heterogeneity along the radial and/or vertical direc- 
tion. For instance, as shown in the present paper, the 
flow of wormlike micelles becomes inhomogeneous first 
along the radial direction r because of the shear-banding 
instability and subsequently along the vertical direction 
z due to secondary flows triggered by an elastic insta- 
bility. It is also possible for some complex fluids to first 
become heterogeneous along z due to a constitutive in- 
stability similar to shear banding and coined "vorticity 
banding"—. On the other hand, yield stress fluids, such 
as emulsions, gels, muds, clay suspensions, or colloidal 
pastes, arc known to present shear localization along r, 
in which a liquidlikc band coexists with a solidlike re- 
gion under simple shear—. In any case, our technique 
is very well suited to investigate these various materials 



in the Taylor- Couette geometry. Ultimately we plan on 
extending our technique to other flows commonly used 
in rheological research, such as the cone-and-plate and 
plate-and-plate geometries. 
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